Identifying Residues for Substrate Recognition in Human GPAT4 by Molecular Dynamics Simulations.

Glycerol-3-phosphate acyltransferase (GPAT) catalyzes the first step in triacylglycerol synthesis. Understanding its substrate recognition mechanism may help to design drugs to regulate the production of glycerol lipids in cells. In this work, we investigate how the native substrate, glycerol-3-phosphate (G3P), and palmitoyl-coenzyme A (CoA) bind to the human GPAT isoform GPAT4 via molecular dynamics simulations (MD). As no experimentally resolved GPAT4 structure is available, the AlphaFold model is employed to construct the GPAT4-substrate complex model. Using another isoform, GPAT1, we demonstrate that once the ligand binding is properly addressed, the AlphaFold complex model can deliver similar results to the experimentally resolved structure in MD simulations. Following the validated protocol of complex construction, we perform MD simulations using the GPAT4-substrate complex. Our simulations reveal that R427 is an important residue in recognizing G3P via a stable salt bridge, but its motion can bring the ligand to different binding hotspots on GPAT4. Such high flexibility can be attributed to the flexible region that exists only on GPAT4 and not on GPAT1. Our study reveals the substrate recognition mechanism of GPAT4 and hence paves the way towards designing GPAT4 inhibitors.

Multiomics profiling of the therapeutic effect of Dan-deng-tong-nao capsule on cerebral ischemia-reperfusion injury.

BACKGROUND: Stroke is a complex physiological process associated with intestinal flora dysbiosis and metabolic disorders. Dan-deng-tong-nao capsule (DDTN) is a traditional Chinese medicine used clinically to treat cerebral ischemia-reperfusion injury (CIRI) for many years. However, little is known about the effects of DDTN in the treatment of CIRI from the perspective of gut microbiota and metabolites. PURPOSE: This study aimed to investigate the regulatory roles of DDTN in endogenous metabolism and gut microbiota in CIRI rats, thus providing a basis for clinical rational drug use and discovering natural products with potential physiological activities in DDTN for the treatment of CIRI. METHODS: The chemical composition of DDTN in vitro and in vivo was investigated using ultra-high performance liquid chromatography-high resolution mass spectrometry (UHPLCHRMS), followed by target prediction using reverse molecular docking. Secondly, a biological evaluation of DDTN ameliorating neural damage in CIRI was performed at the whole animal level. Then, an integrated omics approach based on UHPLCHRMS and 16S rRNA sequencing was proposed to reveal the anti-CIRI effect and possible mechanism of DDTN. Finally, exploring the intrinsic link between changes in metabolite profiles, changes in the intestinal flora, and targets of components to reveal DDTN for the treatment of CIRI. RESULTS: A total of 112 chemical components of DDTN were identified in vitro and 10 absorbed constituents in vivo. The efficacy of DDTN in the treatment of CIRI was confirmed by alleviating cerebral infarction and neurological deficits. After the DDTN intervention, 21 and 26 metabolites were significantly altered in plasma and fecal, respectively. Based on the fecal microbiome, a total of 36 genera were enriched among the different groups. Finally, the results of the network integration analysis showed that the 10 potential active ingredients of DDTN could mediate the differential expression of 24 metabolites and 6 gut microbes by targeting 25 target proteins. CONCLUSION: This study was the first to outline the landscapes of metabolites as well as gut microbiota regulated by DDTN in CIRI rats using multi-omics data, and comprehensively revealed the systematic relationships among ingredients, targets, metabolites, and gut microbiota, thus providing new perspectives on the mechanism of DDTN in the treatment of CIRI.

The 3-hydroxyacyl-CoA dehydratase 1/2 form complex with trans-2-enoyl-CoA reductase involved in substrates transfer in very long chain fatty acid elongation.

Very long-chain fatty acids (VLCFAs) are fatty acids with a carbon chain length greater than 18 carbons (>C18) and exhibit various functions, such as in skin barrier formation, liver homeostasis, myelin maintenance, spermatogenesis, retinal function, and anti-inflammation. VLCFAs are absorbed by dietary or elongated from endogenous hexadecanoyl acids (C16). Similar to long-chain fatty acid synthesis, VLCFAs elongation begins with acyl-CoA and malonyl-CoA as sources, and the length of the acyl chain is extended by two carbon units in each cycle. However, the VLCFAs elongation machinery is located in ER membrane and consists of four components, FA elongase (ELOVL), 3-ketoacyl-CoA reductase (KAR), 3-hydroxyacyl-CoA dehydratase (HACD), and trans-2-enoyl-CoA reductase (TECR), which is different with the long-chain fatty acid machinery fatty acid synthase (FAS) complex. Although the critical components in the elongation cycle are identified, the detailed catalytic and regulation mechanisms are still poorly understood. Here, we focused on the structural and biochemical analysis of TECR-associated VLCFA elongation reactions. Firstly, we identified a stable complex of human HACD2-TECR based on extensive in vitro characterizations. Combining computational modeling and biochemical analysis, we confirmed the critical interactions between TECR and HACD1/2. Then, we proposed the putative substrate binding sites and catalytic residues for TECR and HACD2. Besides, we revealed the structural similarities of HACD with ELOVLs and proposed the possible competition mechanism of TECR-associated complex formation.

Fine-tuning activation specificity of G-protein-coupled receptors via automated path searching.

Physics-based simulation methods can grant atomistic insights into the molecular origin of the function of biomolecules. However, the potential of such approaches has been hindered by their low efficiency, including in the design of selective agonists where simulations of myriad protein-ligand combinations are necessary. Here, we describe an automated input-free path searching protocol that offers (within 14 d using Graphics Processing Unit servers) a minimum free energy path (MFEP) defined in high-dimension configurational space for activating sphingosine-1-phosphate receptors (S1PRs) by arbitrary ligands. The free energy distributions along the MFEP for four distinct ligands and three S1PRs reached a remarkable agreement with Bioluminescence Resonance Energy Transfer (BRET) measurements of G-protein dissociation. In particular, the revealed transition state structures pointed out toward two S1PR3 residues F263/I284, that dictate the preference of existing agonists CBP307 and BAF312 on S1PR1/5. Swapping these residues between S1PR1 and S1PR3 reversed their response to the two agonists in BRET assays. These results inspired us to design improved agonists with both strong polar head and bulky hydrophobic tail for higher selectivity on S1PR1. Through merely three in silico iterations, our tool predicted a unique compound scaffold. BRET assays confirmed that both chiral forms activate S1PR1 at nanomolar concentration, 1 to 2 orders of magnitude less than those for S1PR3/5. Collectively, these results signify the promise of our approach in fine agonist design for G-protein-coupled receptors.

Structural insights into the activation and inhibition of CXC chemokine receptor 3.

The chemotaxis of CD4+ type 1 helper cells and CD8+ cytotoxic lymphocytes, guided by interferon-inducible CXC chemokine 9-11 (CXCL9-11) and CXC chemokine receptor 3 (CXCR3), plays a critical role in type 1 immunity. Here we determined the structures of human CXCR3-DNGi complexes activated by chemokine CXCL11, peptidomimetic agonist PS372424 and biaryl-type agonist VUF11222, and the structure of inactive CXCR3 bound to noncompetitive antagonist SCH546738. Structural analysis revealed that PS372424 shares a similar orthosteric binding pocket to the N terminus of CXCL11, while VUF11222 buries deeper and activates the receptor in a distinct manner. We showed an allosteric binding site between TM5 and TM6, accommodating SCH546738 in the inactive CXCR3. SCH546738 may restrain the receptor at an inactive state by preventing the repacking of TM5 and TM6. By revealing the binding patterns and the pharmacological properties of the four modulators, we present the activation mechanisms of CXCR3 and provide insights for future drug development.

Active DHEA uptake in the prostate gland correlates with aggressive prostate cancer.

Strategies for patient stratification and early intervention are required to improve clinical benefits for patients with prostate cancer. Here, we found that active DHEA utilization in the prostate gland correlated with tumor aggressiveness at early disease stages, and 3ßHSD1 inhibitors were promising for early intervention. [3H]-labeled DHEA consumption was traced in fresh prostatic biopsies ex vivo. Active DHEA utilization was more frequently found in patients with metastatic disease or therapy-resistant disease. Genetic and transcriptomic features associated with the potency of prostatic DHEA utilization were analyzed to generate clinically accessible approaches for patient stratification. UBE3D, by regulating 3ßHSD1 homeostasis, was discovered to be a regulator of patient metabolic heterogeneity. Equilin suppressed DHEA utilization and inhibited tumor growth as a potent 3ßHSD1 antagonist, providing a promising strategy for the early treatment of aggressive prostate cancer. Overall, our findings indicate that patients with active prostatic DHEA utilization might benefit from 3ßHSD1 inhibitors as early intervention.

Approach1es for evolutionary, biochemical, and structural analysis of bacterial steroid 5α-reductases.

Steroid 5α-reductases (SRD5As), also known as 3-oxo-5α-steroid 4-dehydrogenases, are essential membrane-bound enzymes involved in steroid metabolism. Belonging to the NADPH-dependent oxidoreductase family, 5α-reductases catalyze steroids with 3-oxo-Δ4 structure, such as testosterone or progesterone, to produce their corresponding 3-oxo-5α steroids, which are necessary for a variety of physiological and pathological activities. Despite their significance, SRD5A structures are still in short supply to date. Here we describe a protocol for expression, purification, crystallization, structural determination, and functional analysis of PbSRD5A, the 5α-reductase from Proteobacteria bacterium sharing high sequence identity with human SRD5A1 and SRD5A2 isozymes, which we have recently structurally characterized using a lipidic cubic phase approach. Application of similar methods to other 5α-reductase isozymes will lead to breakthroughs in the understanding of the structure, function, and mechanism of oxidoreductases implicated in steroid metabolism.

Mechanism of agonist-induced activation of the human itch receptor MRGPRX1.

Mas-related G-protein-coupled receptors X1-X4 (MRGPRX1-X4) are 4 primate-specific receptors that are recently reported to be responsible for many biological processes, including itch sensation, pain transmission, and inflammatory reactions. MRGPRX1 is the first identified human MRGPR, and its expression is restricted to primary sensory neurons. Due to its dual roles in itch and pain signaling pathways, MRGPRX1 has been regarded as a promising target for itch remission and pain inhibition. Here, we reported a cryo-electron microscopy (cryo-EM) structure of Gq-coupled MRGPRX1 in complex with a synthetic agonist compound 16 in an active conformation at an overall resolution of 3.0 Å via a NanoBiT tethering strategy. Compound 16 is a new pain-relieving compound with high potency and selectivity to MRGPRX1 over other MRGPRXs and opioid receptor. MRGPRX1 was revealed to share common structural features of the Gq-mediated receptor activation mechanism of MRGPRX family members, but the variable residues in orthosteric pocket of MRGPRX1 exhibit the unique agonist recognition pattern, potentially facilitating to design MRGPRX1-specific modulators. Together with receptor activation and itch behavior evaluation assays, our study provides a structural snapshot to modify therapeutic molecules for itch relieving and analgesia targeting MRGPRX1.

An LC-MS-based workflow measures the export activity of S1P transporters.

Sphingosine-1-phosphate (S1P) is an active signaling metabolite synthesized by blood cells, exported into blood stream, and can trigger many downstream signaling pathways with disease implications. Understanding how S1P is transported is of great values for dissecting the function of S1P, but most existing methods for measuring S1P transporter activity use radioactive substrates or involve multiple workup steps, hindering their broader uses. In this study, we develop a workflow combining sensitive LC-MS measurement and a cell-based transporter protein system to measure the export activity of S1P transporter proteins. Our workflow demonstrated good applications in studying different S1P transporters SPNS2 and MFSD2B, WT and mutated protein, and different protein substrates. In summary, we provide a simple yet versatile workflow for measuring the export activity of S1P transporters, facilitating future studies of S1P transport mechanism and drug development.

Structural basis of the substrate recognition and inhibition mechanism of Plasmodium falciparum nucleoside transporter PfENT1.

By lacking de novo purine biosynthesis enzymes, Plasmodium falciparum requires purine nucleoside uptake from host cells. The indispensable nucleoside transporter ENT1 of P. falciparum facilitates nucleoside uptake in the asexual blood stage. Specific inhibitors of PfENT1 prevent the proliferation of P. falciparum at submicromolar concentrations. However, the substrate recognition and inhibitory mechanism of PfENT1 are still elusive. Here, we report cryo-EM structures of PfENT1 in apo, inosine-bound, and inhibitor-bound states. Together with in vitro binding and uptake assays, we identify that inosine is the primary substrate of PfENT1 and that the inosine-binding site is located in the central cavity of PfENT1. The endofacial inhibitor GSK4 occupies the orthosteric site of PfENT1 and explores the allosteric site to block the conformational change of PfENT1. Furthermore, we propose a general "rocker switch" alternating access cycle for ENT transporters. Understanding the substrate recognition and inhibitory mechanisms of PfENT1 will greatly facilitate future efforts in the rational design of antimalarial drugs.

Exploring the Chemical Space of CYP17A1 Inhibitors Using Cheminformatics and Machine Learning.

Cytochrome P450 17A1 (CYP17A1) is one of the key enzymes in steroidogenesis that produces dehydroepiandrosterone (DHEA) from cholesterol. Abnormal DHEA production may lead to the progression of severe diseases, such as prostatic and breast cancers. Thus, CYP17A1 is a druggable target for anti-cancer molecule development. In this study, cheminformatic analyses and quantitative structure-activity relationship (QSAR) modeling were applied on a set of 962 CYP17A1 inhibitors (i.e., consisting of 279 steroidal and 683 nonsteroidal inhibitors) compiled from the ChEMBL database. For steroidal inhibitors, a QSAR classification model built using the PubChem fingerprint along with the extra trees algorithm achieved the best performance, reflected by the accuracy values of 0.933, 0.818, and 0.833 for the training, cross-validation, and test sets, respectively. For nonsteroidal inhibitors, a systematic cheminformatic analysis was applied for exploring the chemical space, Murcko scaffolds, and structure-activity relationships (SARs) for visualizing distributions, patterns, and representative scaffolds for drug discoveries. Furthermore, seven total QSAR classification models were established based on the nonsteroidal scaffolds, and two activity cliff (AC) generators were identified. The best performing model out of these seven was model VIII, which is built upon the PubChem fingerprint along with the random forest algorithm. It achieved a robust accuracy across the training set, the cross-validation set, and the test set, i.e., 0.96, 0.92, and 0.913, respectively. It is anticipated that the results presented herein would be instrumental for further CYP17A1 inhibitor drug discovery efforts.

Structural insights into the covalent regulation of PAPP-A activity by proMBP and STC2.

Originally discovered in the circulation of pregnant women as a protein secreted by placental trophoblasts, the metalloprotease pregnancy-associated plasma protein A (PAPP-A) is also widely expressed by many other tissues. It cleaves insulin-like growth factor-binding proteins (IGFBPs) to increase the bioavailability of IGFs and plays essential roles in multiple growth-promoting processes. While the vast majority of the circulatory PAPP-A in pregnancy is proteolytically inactive due to covalent inhibition by proform of eosinophil major basic protein (proMBP), the activity of PAPP-A can also be covalently inhibited by another less characterized modulator, stanniocalcin-2 (STC2). However, the structural basis of PAPP-A proteolysis and the mechanistic differences between these two modulators are poorly understood. Here we present two cryo-EM structures of endogenous purified PAPP-A in complex with either proMBP or STC2. Both modulators form 2:2 heterotetramer with PAPP-A and establish extensive interactions with multiple domains of PAPP-A that are distal to the catalytic cleft. This exosite-binding property results in a steric hindrance to prevent the binding and cleavage of IGFBPs, while the IGFBP linker region-derived peptides harboring the cleavage sites are no longer sensitive to the modulator treatment. Functional investigation into proMBP-mediated PAPP-A regulation in selective intrauterine growth restriction (sIUGR) pregnancy elucidates that PAPP-A and proMBP collaboratively regulate extravillous trophoblast invasion and the consequent fetal growth. Collectively, our work reveals a novel covalent exosite-competitive inhibition mechanism of PAPP-A and its regulatory effect on placental function.

Prior knowledge facilitates low homologous protein secondary structure prediction with DSM distillation.

MOTIVATION: Protein secondary structure prediction (PSSP) is one of the fundamental and challenging problems in the field of computational biology. Accurate PSSP relies on sufficient homologous protein sequences to build the multiple sequence alignment (MSA). Unfortunately, many proteins lack homologous sequences, which results in the low quality of MSA and poor performance. In this article, we propose the novel dynamic scoring matrix (DSM)-Distil to tackle this issue, which takes advantage of the pretrained BERT and exploits the knowledge distillation on the newly designed DSM features. Specifically, we propose the DSM to replace the widely used profile and PSSM (position-specific scoring matrix) features. DSM could automatically dig for the suitable feature for each residue, based on the original profile. Namely, DSM-Distil not only could adapt to the low homologous proteins but also is compatible with high homologous ones. Thanks to the dynamic property, DSM could adapt to the input data much better and achieve higher performance. Moreover, to compensate for low-quality MSA, we propose to generate the pseudo-DSM from a pretrained BERT model and aggregate it with the original DSM by adaptive residue-wise fusion, which helps to build richer and more complete input features. In addition, we propose to supervise the learning of low-quality DSM features using high-quality ones. To achieve this, a novel teacher-student model is designed to distill the knowledge from proteins with high homologous sequences to that of low ones. Combining all the proposed methods, our model achieves the new state-of-the-art performance for low homologous proteins. RESULTS: Compared with the previous state-of-the-art method 'Bagging', DSM-Distil achieves an improvement about 5% and 7.3% improvement for proteins with MSA count ≤30 and extremely low homologous cases, respectively. We also compare DSM-Distil with Alphafold2 which is a state-of-the-art framework for protein structure prediction. DSM-Distil outperforms Alphafold2 by 4.1% on extremely low-quality MSA on 8-state secondary structure prediction. Moreover, we release a large-scale up-to-date test dataset BC40 for low-quality MSA structure prediction evaluation. AVAILABILITY AND IMPLEMENTATION: BC40 dataset: https://drive.google.com/drive/folders/15vwRoOjAkhhwfjDk6-YoKGf4JzZXIMC. HardCase dataset: https://drive.google.com/drive/folders/1BvduOr2b7cObUHy6GuEWk-aUkKJgzTUv. Code: https://github.com/qinwang-ai/DSM-Distil.

Inhibiting 3ßHSD1 to eliminate the oncogenic effects of progesterone in prostate cancer.

Prostate cancer continuously progresses following deprivation of circulating androgens originating from the testis and adrenal glands, indicating the existence of oncometabolites beyond androgens. In this study, mass-spectrometry-based screening of clinical specimens and a retrospective analysis on the clinical data of prostate cancer patients indicate the potential oncogenic effects of progesterone in patients. High doses of progesterone activate canonical and non-canonical androgen receptor (AR) target genes. Physiological levels of progesterone facilitate cell proliferation via GATA2. Inhibitors of 3ß-hydroxysteroid dehydrogenase 1 (3ßHSD1) has been discovered and shown to suppress the generation of progesterone, eliminating its transient and accumulating oncogenic effects. An increase in progesterone is associated with poor clinical outcomes in patients and may be used as a predictive biomarker. Overall, we demonstrate that progesterone acts as an oncogenic hormone in prostate cancer, and strategies to eliminate its oncogenic effects may benefit prostate cancer patients.

Structural insights into sphingosine-1-phosphate receptor activation.

As a critical sphingolipid metabolite, sphingosine-1-phosphate (S1P) plays an essential role in immune and vascular systems. There are five S1P receptors, designated as S1PR1 to S1PR5, encoded in the human genome, and their activities are governed by endogenous S1P, lipid-like S1P mimics, or nonlipid-like therapeutic molecules. Among S1PRs, S1PR1 stands out due to its nonredundant functions, such as the egress of T and B cells from the thymus and secondary lymphoid tissues, making it a potential therapeutic target. However, the structural basis of S1PR1 activation and regulation by various agonists remains unclear. Here, we report four atomic resolution cryo-electron microscopy (cryo-EM) structures of Gi-coupled human S1PR1 complexes: bound to endogenous agonist d18:1 S1P, benchmark lipid-like S1P mimic phosphorylated Fingolimod [(S)-FTY720-P], or nonlipid-like therapeutic molecule CBP-307 in two binding modes. Our results revealed the similarities and differences of activation of S1PR1 through distinct ligands binding to the amphiphilic orthosteric pocket. We also proposed a two-step "shallow to deep" transition process of CBP-307 for S1PR1 activation. Both binding modes of CBP-307 could activate S1PR1, but from shallow to deep transition may trigger the rotation of the N-terminal helix of Gαi and further stabilize the complex by increasing the Gαi interaction with the cell membrane. We combine with extensive biochemical analysis and molecular dynamic simulations to suggest key steps of S1P binding and receptor activation. The above results decipher the common feature of the S1PR1 agonist recognition and activation mechanism and will firmly promote the development of therapeutics targeting S1PRs.

Crystal structure of steroid reductase SRD5A reveals conserved steroid reduction mechanism.

Steroid hormones are essential in stress response, immune system regulation, and reproduction in mammals. Steroids with 3-oxo-Δ4 structure, such as testosterone or progesterone, are catalyzed by steroid 5α-reductases (SRD5As) to generate their corresponding 3-oxo-5α steroids, which are essential for multiple physiological and pathological processes. SRD5A2 is already a target of clinically relevant drugs. However, the detailed mechanism of SRD5A-mediated reduction remains elusive. Here we report the crystal structure of PbSRD5A from Proteobacteria bacterium, a homolog of both SRD5A1 and SRD5A2, in complex with the cofactor NADPH at 2.0 Å resolution. PbSRD5A exists as a monomer comprised of seven transmembrane segments (TMs). The TM1-4 enclose a hydrophobic substrate binding cavity, whereas TM5-7 coordinate cofactor NADPH through extensive hydrogen bonds network. Homology-based structural models of HsSRD5A1 and -2, together with biochemical characterization, define the substrate binding pocket of SRD5As, explain the properties of disease-related mutants and provide an important framework for further understanding of the mechanism of NADPH mediated steroids 3-oxo-Δ4 reduction. Based on these analyses, the design of therapeutic molecules targeting SRD5As with improved specificity and therapeutic efficacy would be possible.

A Glimpse of the Structural Biology of the Metabolism of Sphingosine-1-Phosphate.

As a key sphingolipid metabolite, sphingosine-1-phosphate (S1P) plays crucial roles in vascular and immune systems. It regulates angiogenesis, vascular integrity and homeostasis, allergic responses, and lymphocyte trafficking. S1P is interconverted with sphingosine, which is also derived from the deacylation of ceramide. S1P levels and the ratio to ceramide in cells are tightly regulated by its metabolic pathways. Abnormal S1P production causes the occurrence and progression of numerous severe diseases, such as metabolic syndrome, cancers, autoimmune disorders such as multiple sclerosis, and kidney and cardiovascular diseases. In recent years, huge advances on the structure of S1P metabolic pathways have been accomplished. In this review, we have got a glimpse of S1P metabolism through structural and biochemical studies of: sphingosine kinases, S1P transporters and S1P receptors, and the development of therapeutics targeting S1P signaling. The progress we summarize here could provide fresh perspectives to further the exploration of S1P functions and facilitate the development of therapeutic molecules targeting S1P signaling with improved specificity and therapeutic effects.